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Abstract 

Based on the standard cosmological model, we calculate the correction to the rate of 
two-photon 2s <-» Is transitions in the hydrogen atom under primordial hydrogen plasma 
recombination conditions that arises when the induced transitions under equilibrium back- 
ground radiation with a blackbody spectrum and plasma recombination radiation are 
taken into account. 
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1 Introduction 



The cosmological plasma recombination plays an important role in producing the cosmic 
microwave background radiation (CMBR) anisotropy, since it determines the time when 
the optical depth of the Universe for Compton scattering becomes less than unity. The 
recombination model suggested by Zel'dovich et al. (1968) and Peebles (1968) for the 
hydrogen plasma and generalized by Seager et al. (1999, 2000) to the hydrogen-helium 
plasma allows an accuracy of at least 1% to be achieved when calculating the cosmological 
recombination. Therefore, in the last few years, the efforts of researchers in this field have 
been directed at studying the effects that can affect the primordial plasma recombination 
kinetics at a level of up to several percent (Leung et al. 2004; Dubrovich and Grachev 
2005; Chluba and Sunyaev 2006; Novosyadlyj 2006). 

The recombination kinetics depends mainly on two processes: the two-photon 2s <-> 
Is transition and the Ly-a photon escape from the 2p — > Is line profile due to the 
cosmological redshift (Zel'dovich et al. 1968; Peebles 1968; Seager et al. 1999, 2000). 

In this paper, we discuss the correction to the rate of two-photon 2s <-> Is transi- 
tions due to the effect of induced transitions. Such an effect was first discussed by Chluba 
and Sunyaev (2006), who showed that including it leads to the relative change in the 
ionization fraction by one percent for epochs z = 900 — 1200. This is important, since, in 
turn, it can lead to changes in the CMBR temperature fluctuation spectrum, which can 
be measured in planned experiments (Planck and others). 

In this paper, special emphasis was placed on studying the Is — > 2s transi- 
tions (/iz/ 2 i ~ 10.2 eV) related to two-photon absorption. In addition to the study by 
Chluba and Sunyaev (2006), we took into account the fact that the nonequilibrium Ly-a 
photons produced under primordial hydrogen plasma recombination are shifted to the 
lower frequencies. For the recombination period, the latter are limited by the frequency 
v r 2 f = ^21 (1 + z end )/(l + z begin ) ~ 0.6i/ 2 i , where z begin ~ 1600 and z end ~ 900 define the 
period in which more than 95% of the recombination Ly-a photons were produced. As 
a result of this process, the occupation number of the primordial plasma recombination 
radiation photons (reddened Ly-a photons) at the recombination epoch in the frequency 
range [O.6Z/21; ^21] exceeds the occupation number of the equilibrium background photons 
by 1-4 orders of magnitude (see Fig. 1). This, in turn, results in a significant difference 
between the Is — > 2s reaction rate and that from Chluba and Sunyaev (2006), who as- 
sumed the presence of only equilibrium background radiation with a blackbody spectrum 
in the above frequency range. 

2 Uncompensated 2s <-> Is transitions 

Calculating the correction to the rate of two-photon 2s <-> Is transitions is based on 
calculating the uncompensated electron flow J2 S u between the 2s and Is levels. The 
uncompensated electron flow J 2s i s is described by the expression 

J2sls — P2sls%2s — Pls2s%ls (l) 

Here, P2sis [s -1 ] is the rate of two-photon 2s — > Is transitions including the spontaneous 
and induced transitions, Pi S 2 S is the rate of two-photon Is — > 2s transitions including 
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the induced transitions, X2 S and x\ s are the relative populations of the 2s and Is states 
defined by the expressions 



x 2s = N 2s /N H , x u = N ls /N H (2) 

where N 2s is the number density of 2s-state atoms, N hs is the number density of ground- 
state atoms, and Nh is the total number density of neutral and ionized hydrogen atoms. 

Let us define the rate of uncompensated 2s — * Is transitions A% sls (the superscript 
u stands for uncompensated) as 

^2sls = J2sls/X2s = P2.sU — P\s2s , (3) 

X2s 

A2 S i s can be comparable to the probability A 2s i s — 8.22 s -1 of spontaneous two-photon 
2s — > Is decays that determines the rate of 2s — > Is transitions in the absence of an 
external radiation field (i.e., in a vacuum). In turn, the probability A 2s i s of spontaneous 
2s — > Is transitions in the hydrogen atom is defined by the expression (Spitzer and 
Greenstein 1951; Zon and Rapoport 1968) 

A2su = y f <Kv)dv ( 4 ) 

where A = 9a 6 cRy/2 10 ~ 4.3663 s" 1 is the dimensional normalization constant, a is the 
fine-structure constant, Ry is the Rydberg constant for hydrogen, c is the speed of light, 
and the function (f>(v/v2i) describes the spontaneous two-photon decay spectrum (Spitzer 
and Greenstein 1951; Zon and Rapoport 1968). The value of <f>(y)dy is proportional to 
the spontaneous emission probability of a pair of photons one of which has frequency v in 
the range [i^iy; 1*212/ + V2\dy\ and the other has frequency v' = (i/ 2 i — v). The factor 1/2 
arises, because each pair of photons is taken into account twice when integrating in (4). 

Calculating the rate of uncompensated transitions AV^ sXs requires knowledge of P2 S i s 
and P\ S 2s- The transition rates P2 S i s and P\ S 2s depend on the parameters of the radiation 
field with which the primordial plasma interacts. This dependence is taken into account 
by the expressions (Rapoport et al. 1978) 

p -« = tW£> (i+,m)(i+,( M£)' <5) 

P - = tH£)^w(£), (6) 

where r){y) is the occupation number at frequency v in the radiation field with which the 
primordial plasma interacts. 

Thus, taking Eqs. (3)- (6) into consideration, we can represent A\ ls in the following 
form containing an explicit dependence on the occupation numbers r\: 



;i + ^))(i + ^'))-^M^)— 

X2s 



Basically, Eq. (7) is the collision integral for photons and hydrogen atoms. 
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To simplify the subsequent calculations, let us pass to the following expression 
for the rate of uncompensated transitions with integration over half the range of the 
2s — > Is transition, which is valid, since the integrands are symmetric with respect to the 
substitution v <-> (v 21 — v): 



AU 

/1 2sls 



An 



1/2 



V 



1 + 77(1/)) (1 + T){v')) - T){v)t){v') — 

%2s 



^21 



(8) 



Although the integration in Eq. (8) is over the frequency range [0; 1*21/2], the 
argument v' = {v 2 \ — v) varies over the range [u 2 i /2; 1/21], i.e., the integrand depends on 
the occupation numbers f] in the full transition frequency range [0; v 2 i\- 

Since the frequency of more than 95% of the Ly-a photons does not decrease to 
1*21/2 in the recombination period and since the perturbations of the occupation numbers 
related to the transitions between excited states of the hydrogen atom (transition frequen- 
cies < 1*21/2) are no larger than 10~ 5 (see, e.g., Dubrovich and Grachev 2004; Kholupenko 
et al. 2005), the occupation number rj(u) in the frequency range < v < 1*21/2 is described 
with a relative accuracy of at least 10~ 3 by the Planck distribution 



cxp 



hv 
kTf 



- 1 



(9) 



where ks is the Boltzmann constant, the parameter T = T (l + z) is the CMBR tem- 
perature at epoch z, and T = 2.726 K is the CMBR temperature at the present epoch. 
Accordingly, the following approximate equality is valid for A 2sls : 



AU 

^2sls 



Ar 



1/2 



V 
"21 



Xl s 
X2s 



V 



\V2\ 



(10) 



The second approximation that allows Eq. (8) to be simplified stems from the fact that 
the occupation numbers r\{y r ) are no larger in absolute value than 10 -5 (i.e., r\{y') <^ 1) 
in the frequency range 1*21/2 < v' < v 2 \ at the recombination epoch. Therefore, r]{y') in 
the first term in the integrand in (10) may be ignored compared to unity. At the same 
time, r]{v') cannot be ignored in the second term, since it appears in Eq. (10) in the 
product with a large factor X\ s /x 2s 3> 1. Taking this into account, we obtain a simplified 
expression for A2 Sls : 



A U - An 



1/2 



V 



X2s 



d 



V 



V21 



Let us introduce the function 



Rly) = 1 - 



r](v')rf(v) x ls 

1 + X 2t 



(11) 



(12) 



that characterizes the 2s <-> Is reaction rate at certain emitted photon frequencies v and 
v' = (i*2i — v). If the rate of the reverse reaction is negligible compared to the rate of the 
direct reaction (e.g., for late epochs z < 1000), then R ~ 1. If, alternatively, the direct 
and reverse reactions compensate each other (e.g., at early epochs z > 1500 when the 
populations are close to the equilibrium ones), then R ~ 0. 
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Using definition (12), we can write Eq. (11) in a compact form: 



r*fe) (i+ " v))flMd fe) <i3) 



The optical depth in the Ly-a line at the recombination epoch is 10 5 — 10 8 , which 
leads to the following relation between the populations x is , x 2p and the occupation number 
in the Lya line: ^(^21) = x 2p /3xi s . Since the condition x 2p = 3x 2s is satisfied with a 
relative accuracy of at least 10~ 2 for epochs z > 800 (Grachev and Dubrovich 1991), we 
obtain the following expression for 77(^21) (see also Peebles 1968): 

= — (14) 
Substituting Eqs. (9) and (14) in Eq. (12) yields 

w-i 1 - <15) 

The occupation number ^(^21) in the Ly-a line can be calculated using the recom- 
bination equation (Peebles 1968). In our notation, this equation is 



x P = -a c N H x + (3 c r](u 21 )(l - x p ) (16) 



where x p = N p /Nn is the hydrogen plasma ionization fraction, a c is the total recombina- 
tion coefficient to all excited states of the hydrogen atom, and j3 c is the total ionization 
coefficient from all excited states. These coefficients are related by 

p c = a c g e (T)ex P (^-^j (17) 

where g e {T) is the partition function of free electrons and E 2 is the ionization energy of 
the hydrogen atom in the 2s state. 

Using (16) and (17), we can write 

,7(1*0 = 7(^°W(1 + *(*)), S(z) = (18) 

where 7(2) is the function that characterizes the plasma deviation from ionization equi- 
librium, 

4 N H ( E 1 



where Ei is the ionization energy of the hydrogen atom in the ground state. Calculating 
the function 7(2) requires prespecifying the dependence x p (z), which is determined by 
solving the problem of primordial plasma recombination through numerical calculation 
(e.g., using recfast; Seager et al. 1999). The calculation of the recombination kinetics 
itself may be performed without including the effect considered here, i.e., in terms of the 
standard recombination model (Peebles 1968; Seager et al. 1999, 2000), since the expected 
relative change in the ionization fraction after applying the correction under study to the 
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rate of two-photon 2s <-> Is decay is ~ 10~ 2 (Chluba and Sunyaev 2006). The function 
7(2) is presented in Fig. 2. In this paper, we will consider 7(2) as a known quantity in 
terms of which all of the soughtfor functions will be expressed. 

Since the value of the function \5(z)\ for epochs z > 1100 does not exceed 2 • 10~ 2 , 
below we will use an approximate equality to determine 77(1*21) (disregarding S(z) compared 
to unity): 

?7(^ 21 ,^) ~ 7(^)77°(zy 21 ,^) (20) 

The last function to be determined to calculate R{y) is the total occupation number 
r]{v') of CMBR photons and recombination Ly-a photons whose frequency decreased from 
the initial value of 1*21 to v' as 

, 1 + z , . 

✓ = , 21TT7 (21) 

where z' is the epoch when the photon frequency was z/ 2 i. Taking into account the 
law of frequency change (21), we obtain 77(1/, 2) = 77(1*21,2') and rf{v',z) = rf{v 2 \,z'). 
Combining these equalities and equality (20) leads to the following formulae to determine 

Vi ' } r 1 {v\z)=r 1 {v^z')=r 1 \v 2l ,z l ) 1 {z')=r l \v\z) 1 {z') (22) 

where the parameter z' depends on v' and z and is defined by the inequality that follows 
from (21): 

*V,*) = (^r)(i + *)-i (23) 

Given (20) and (22), the formula for R{y) can be written as 

A further simplification of this formula is related to the fact that the following 
equality is valid for v' > u 2 i/2 at the recombination epoch: 

77V) / hu \ 

ex P ( ~17~^ ) - 1 (25) 



77 M 'V k BT 
since the approximate equality 

rf{v') ~ exp (-^) , (26) 

holds with a relative accuracy of at least 10~ 5 in this period. Taking this into account, 
we obtain the final formula for R(v): 

R{y) = (1 - l{z')h{z)) (27) 

Substituting this expression in Eq. (13) yields the final expression for the rate of uncom- 
pensated decays 

ALu = A. r (—) (1 + V») (1 - l(z')h(z)) d (—) (28) 



If equilibrium radiation is assumed to be present in the frequency range [0; ^21], then 
the principle of detailed balancing can be used (Chluba and Sunyaev 2006). The rate of 
uncompensated transitions calculated under this assumption is given by the expression 
(in addition, highlighted by the superscript CS): 

A 2 fJ = A (^) (1 + V°(", z)) (1 - l/7(*)) d [^j (29) 

In the standard recombination calculation (Peebles 1968; Seager et al. 1999, 2000), 
no induced transitions are considered, i.e., the occupation number rf(h>) in Eq. (29) is 
disregarded. Thus, the rate of uncompensated transitions can be calculated using the 
formula (the superscript st stands for standard) 

Since the key aspect of this study is allowance for the recombination Ly-a photons 
(because the Is — > 2s reaction rate depends on their number), we calculated the ratio of 
the intensity of the "reddened" recombination radiation to the intensity of the equilibrium 
background radiation with temperature T that must be observed at the current epoch 
(Fig. 3). This quantity is related to the function 7(2) by 

I{Z = ^distortion = V(v) ~ W) 
I(Z = 0)pianck rf{v) 

where z v = v 2 \ jv — 1 is the redshift of the Ly-a photon produced at the recombination 
epoch. Previously, this quantity was calculated by Wong et al. (2005). Our calculations 
show satisfactory agreement with the calculations by Wong et al. (2005). Our calcula- 
tions also agree with those of Grachev and Dubrovich (1991) and Boschan and Biltzinger 
(1998), who calculated the distortion of the CMBR intensity related to the recombination 
radiation in the Ly-a line. 



(31) 



3 Results 

The calculations were performed in terms of the standard ACDM cosmological model with 
the following parameters at the present epoch: the Hubble constant H = 70 km/s/Mpc; 
the relative vacuum-like energy density (cosmological constant) ft\ ~ 0.73; the relative 
nonrelativistic matter density fl m ~ 0.27; the relative baryonic density ~ 0.036; and 
the relative relativistic matter density Q re i ~ 10~ 4 (photons + neutrinos); 

Figure 4 reproduces the results by Chluba and Sunyaev (2006): the spectrum (solid 
curve) of the spontaneous two-photon 2s — > Is decay 0(z//V 2 i); and the total spectrum 
(dotted curve) of the spontaneous and induced two-photon 2s — > Is decays (Chluba and 
Sunyaev 2006), which is symmetric with respect to v 2 i/2 and is defined by the following 
function in the range [0; u 2 i/2]: 

<P md {— )=<!>{— ) (l+V»), (32) 

\ u 21 / \ V 2 \ J 
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calculated for the epoch z = 1500 (the superscript ind stands for induced). 
Figure 5 presents the following calculated spectra: 

(1) The spectrum of the two-photon 2s <-> Is transition including only the spon- 
taneous decays, with the reverse reaction channel being taken into account using the 
principle of detailed balancing: 



This spectrum was calculated under the assumptions made in the standard recombination 
model (Peebles 1968; Seager et al. 1999). In Fig. 5, it is represented by the solid curve. 

(2) The spectrum of the two-photon 2s <-> Is transition including the spontaneous 
and induced transitions, the direct and reverse reaction channels, and under the assump- 
tion that the radiation is equilibrium one in the entire 2s — > Is transition frequency range: 
(p cs [yjvix). This spectrum is also symmetric with respect to u 2 \/2 and is defined by the 
following function in the range [0; 1*21/2]: 



This spectrum was calculated under the assumptions made by Chluba and Sunyaev (2006). 
In Fig. 5, it is represented by the dotted curve. (3) The spectrum of the two-photon 
2s <-> Is transition including the spontaneous and induced transitions, the direct and 
reverse reaction channels, and the nonequilibrium radiation in the Ly-a line and the 
range 0. 61*21 < v' < u 2 \. This spectrum is also symmetric with respect to 1*21/2 and is 
defined by the following function in the range [0; 1/21/2]: 



This spectrum was calculated under the assumptions made here. In Fig. 5, it is repre- 
sented by the dashed curve. 

In contrast to the frequency-integrated reaction rates A^'/^, A^' sls , and A% sls the 
functions s *, <p cs , and 0" presented in Fig. 5 characterize the behavior of the 2s <-> Is 
reaction rate not only as a function of the epoch (z = 1500 - Fig. 5, top; z = 1100 - Fig. 
5, bottom), but also as a function of the emitted photon frequencies v and v 1 = (1*21 — v). 
The latter is important, since the key aspect of this paper is allowance for the effect of 
reddened recombination Ly-a photons, whose number depends significantly on frequency, 
on the 2s <-> Is reaction rate. 

Note also that the spectra (f> st , <p cs , and M shown in Fig. 5 are of interest in their 
own right when calculating the distortion of the CMBR by the recombination photons 
produced during two-photon 2s <-> Is transitions (Zel'dovich et al. 1968; Wong et al. 
2005), since these functions describe the behavior of the source of recombination photons 
for each specific recombination model: (f> st for the standard recombination model (Peebles 
1968; Seager et al. 1999); (f> cs for the modification of the standard model suggested by 
Chluba and Sunyaev (2006); and 0" for the modification of the standard model considered 
here. 




(33) 




(34) 




(35) 
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The fact that the spectrum n , in contrast to the spectrum (f> cs , becomes zero at 
the boundary of the transition frequency range (at v — 0) is related to condition (14). 
The latter leads to full compensation of the induced two-photon 2s — > Is transitions 
by the induced reverse two-photon Is — > 2s transitions at the frequency v' = v 2 \ that 
corresponds to the Ly-a photon absorption. 

Note also that the behavior of the spectra as a function of the epoch z under 
consideration is peculiar in that for epochs z > 1375 (the model-dependent value obtained 
by numerical calculation), the values of <p u {u) are higher than those of c/> s *(i/) (see Fig. 5, 
top), i.e., A2 Sls > A^gig, while for epochs z < 1375, the values of (f> u {v) are lower than 
those of 4> st {v) (see Fig. 5, bottom), i.e., ^2sis < A 2sis- This is because the 2s level for 
z < 1375 is underpopulated with respect to the Is level, given the occupation number 
rj(u'), i.e., 

x 2s /x ls = vfoi) < r]{v') (36) 

This, in turn, is attributable to the nonmonotonic behavior of rj(u') with frequency and 
epoch z under consideration (see Fig. 6). 

Figure 7 shows the relative changes SA CS and 8 A in the rate of uncompensated 
2s <-> Is transitions compared to the standard calculation (Peebles 1968; Seager et al. 
1999) related to the inclusion of induced transitions. The value of SA CS was calculated 
under the assumptions made by Chluba and Sunyaev (2006); the value of 6 A was calcu- 
lated under the assumptions made here. By definition (3), 5A CS and 5 A also describe the 
relative changes in the uncompensated electron flow J 2s is between the 2s and Is levels 
compared to the standard calculation (Peebles 1968; Seager et al. 1999): 



AU,Ob _ A u,st jCS _ jst 

X ACS — ^sls ^sls _ J 2sls J 2sls /of-A 

^sls J 2sls 



(dotted curve) and 



AU _ AU,St J jst 

r A /1 2sls ^sls J 2sls ^2sls /oo\ 

dA ~ J^Tt - JTt 

(dashed curve). 

Figure 8 presents the results of our calculations of the following quantities: 

(1) Figure 8 (top) shows the calculated total flow J r [cm~ 3 s -1 ] of recombining 
primordial hydrogen plasma electrons: 

J r = H(z)(l + z)N H (z)^, (39) 

where H(z) is the Hubble constant at epoch z. 

(2) Figure 8 (bottom) shows the relative change 5J r in the flow of recombining 
electrons compared to the standard recombination model (Peebles 1968; Seager et al. 
1999) related to the inclusion of induced transitions: 

7 — I st 

Mr = (40) 

Figure 9 shows the calculated primordial hydrogen plasma ionization fraction x p 
as a function of redshift (Fig. 9, top) and the relative change in this quantity 5x p = 
(x p — x^) I 'x s p l compared to the standard recombination model (Peebles 1968; Seager et al. 
1999) related to the inclusion of induced transitions (Fig. 9, bottom). 
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4 Conclusions 



The calculations show that when the radiation produced during the primordial hydrogen 
plasma recombination and the cosmological redshift is taken into account, the correction 
5A to the rate of uncompensated 2s <->• Is transitions due to the induced two-photon 
2s <-> Is transitions is a factor of 2 smaller than that calculated under the assumptions 
made by Chluba and Sunyaev (2006) (see Fig. 7) and even has the negative sign at epochs 
z < 1375, which generally leads to a delay of the primordial plasma recombination (see 
Fig. 9) compared with the standard calculation (Peebles 1968; Seager et al. 1999). The 
relative change 5x p in the ionization fraction for epochs z < 1200 is 0.5%. This, in turn, 
leads to the surface of the last scattering being located at a lower z than that obtained in 
the standard calculation of the primordial plasma recombination kinetics (Peebles 1968; 
Seager et al. 1999). Further calculations are required to accurately calculate the location 
of the surface of the last scattering and the CMBR temperature fluctuation spectrum by 
taking into account the modification of the recombination model suggested here. 
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Figure 1: Fig. 1. Total occupation number 77 for the primordial hydrogen plasma re- 
combination radiation and the background equilibrium radiation vs. frequency v (solid 
curves); occupation number 77° for the background equilibrium radiation vs. frequency v 
(dashed curves) . The excess of the occupation numbers 77 above their equilibrium values 
of 77 (i.e., the difference (77 — 77 )) is determined by the number of "reddened" recom- 
bination Ly-a photons. The curves from the top downward correspond to the epochs 
z = 1400, 1200, 1000. 19 




z 



Figure 2: Fig. 2. Function 'y(z) characterizing the plasma deviation from ionization 
equilibrium calculated in terms of the standard ACDM cosmological model. 



13 




100 120 140 160 180 200 



\>bs' ^ m 



Figure 3: Fig. 3. Ratio of the intensity of the primordial hydrogen plasma recombination 
radiation to the intensity of the CMBR (T = 2.726 K) at the current epoch vs. observed 
wavelength X b s - 
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Figure 4: Fig. 4. Spectrum of the spontaneous two-photon decays (solid curve) and 
spectrum of the two-photon decays including the spontaneous and induced decays <p md 
(dotted curve) for the epoch z=1500 (Chluba and Sunyaev 2006). 
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Figure 5: Fig. 5. Spectra of the two-photon 2s <-> Is transitions at various epochs: 
the spectrum (p st calculated under the assumptions of the standard recombination model 
(Peebles 1968; Seager et al. 1999) (solid curve); the spectrum <p cs calculated under the 
assumptions made by Chluba and Sunyaev (2006) (dotted curve); and the spectrum (f) u 
calculated under the assumptions made here (dashed curve). 
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Figure 6: Fig. 6. Dynamics of the recombination photon occupation number 77: the curves 
for epochs z = 1000 — 1500 are presented. 



17 



0,05 



< 
to 



0,04 



0,03 



0,02 - 



0,01 







-0,01 - 



-0,02 



800 




1600 



Figure 7: Fig. 7. Relative change in the rate of uncompensated 2s <-> Is transitions com- 
pared to the rate of uncompensated 2s <-> Is transitions calculated under the assumptions 
made in the standard model (Peebles 1968; Seager et al. 1999): 5A CS is the result ob- 
tained under the assumptions made by Chluba and Sunyaev (2006) (dotted curve); and 
5 A is the result obtained under the assumptions made here (dashed curve). 
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Figure 8: Fig. 8. Top: Flow of irretrievably recombining primordial hydrogen plasma 
electrons J r . Bottom: The relative change in the flow of recombining electrons 5J r com- 
pared to the standard recombination model (Peebles 1968; Seager et al. 1999) related to 
the inclusion of induced transitions. 
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Figure 9: Fig. 9. Top: Primordial hydrogen plasma ionization fraction x p vs. redshift 
z. Bottom: The relative change 5x p = (x p — xp*)/Xp in the ionization fraction compared 
to the standard recombination model (Peebles 1968; Seager et al. 1999) related to the 
inclusion of induced transitions. 
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